We used matched filter detection and multiple-event relocation techniques to characterize the spatiotemporal evolution of the sequence. Our analysis is from the 14 closest seismic stations to the earthquake sequence, which included seven permanent stations from the Montana Regional Seismic Network, one permanent station from the ANSS backbone network and three temporary seismic stations deployed by the USGS within four days after the mainshock. A catalog of 685 well-located earthquakes larger than M 1 occurring Between 5 July and 15 October 2017 were relocated using a hypocentroid decomposition (HD) multiple-event relocation approach. The resulting dataset had an average epicentral and depth uncertainties (90% confidence) on the order of 1 km. Using match filtering of the station waveforms for each of these relocated events we were able to detect four foreshocks in the three days prior to, and 3005 aftershocks in the three weeks following the mainshock.